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The vortex streets produced by a napping foil of span-to-chord aspect ratio of 4:1 are 
studied in a hydrodynamic tunnel experiment. In particular, the mechanisms giving rise 
to the symmetry breaking of the reverse Benard-von Karman vortex street that char- 
acterizes fish-like swimming and forward flapping flight are examined. Two-dimensional 
particle image velocimetry measurements in the mid-plane perpendicular to the span axis 
of the foil are used to characterize the different flow regimes. The deflection angle of the 
mean jet flow with respect to the horizontal observed in the average velocity field is used 
as a measure of the asymmetry of the vortex street. Time series of the vorticity field are 
used to calculate the advection velocity of the vortices with respect to the free-stream, 
defined as the phase velocity [/phase as well as the circulation Y of each vortex and the 
spacing £ between consecutive vortices in the near wake. The observation that the sym- 
metry breaking results from the formation of a dipolar structure from each couple of 
counter-rotating vortices shed on each flapping period serves as starting point to build a 
model for the symmetry breaking threshold. A symmetry breaking criterion based on the 
relation between the phase velocity of the vortex street and an idealized self-advection 
velocity of two consecutive counter-rotating vortices in the near wake is established. The 
predicted threshold for symmetry breaking accounts well for the deflected wake regimes 
observed in the present experiments and may be useful to explain other experimental 
and numerical observations of similar deflected propulsive vortex streets reported in the 
literature. 



1. Introduction 

Flapping-based propulsive systems, either natural or man-made, are often discussed in 
terms of the Strouhal number, defined as the product of the flap ping frequency f and 
amplitude A divid ed by the cruising speed Uq, i.e. St a — fA/Uo ijAnderson et al. 19981 : 



Taylor et a^.l 12003^ . Another crucial parameter in these problems is the aspect ratio of 



the flapping body, because it determines to what extent a quasi-two-dimensional (Q2D) 
view can capture the main elements needed for an adequate description of the real three- 
dimensional (3D) flow. In particular, in the case of a flapping body propelling itself in 
forward motion, at least two qualitatively different situations have been evidenced from 
flapping foil experiments and numerical simulations: high span-to-chord ratio foils pro- 
duce a propulsive vortex street, the rever se Benard-von Karman (BvK) wake (see e.g. 
lKoochesfahanil[l989T : Anderson et al. 1993), where the most intense vortices are aligned 
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Figure 1. Transitions in the wake of a flapping foil in the Ad vs. St map for Re = 255 (from 
iGodov-Diana et al\\200§i ). Experimental points are labeled as □: BvK wake; ■: aligned vortices 
(2S wake); +: reverse BvK wake; A: deflected reverse BvK street resulting in an asymmetric 
wake. Solid line: transition between BvK and reverse BvK. Dashed line: transition between 
reverse BvK and the asymmetric regime. Typical vorticity fields are shown as inserts on each 
region. 



with the foil span and turn in opposite senses with respect to the natural BvK vortices 
behind a 2D cylinder. A Q2D analysis accounts for the key dynamical features in this 
case where the mean flow has the form of a jet and results in a net propulsive force. As 
the span-to-chord ratio decreases towards unity, 3D effects come into play and modify 
dramatically the structure of the wake. In this case a series of vortex loops (or horse-shoe 
vor tices) are engendered from the vorticity shed from all sid es of the flapping foil (see 
e.g. von Ellenrieder et a?.ll2003l : Buchholz fe Smitsl^OOdEooih . The experiments reported 



here were performed with a 4:1 aspect ratio foil, which is high enough to produce Q2D 
regimes in the near wake. A two-parameter description that permits to vary indepen- 
dently the frequency and amplitude of the oscillatory motion has been shown recently 
I Godov-Diana et al. 20081 ) to be the optimum framework to fully characterize the quasi- 



two-dimensional regimes observed in the wake of a pitching foil. The transition from a 
BvK vortex street to the reverse BvK street characteristic of propulsive regim e s, and the 
symmetry breaking of the reverse BvK street reported in lGodoy-Diana et al. (2008) are 



summarized in the (St, Ad) phase space shown in figure [TJ The Strouhal number and a 
dimensionless amplitude have been defined using a fixed length scale (the foil width D) 
as St = fD/Ug and Ad = A/D, respectively Q. 

In this paper we focus on the reverse BvK regime, attempting to shed some light on 
the physical mechanisms that determine its symmetry breaking (the dashed line in figure 
[U . Asymmetric vor tex streets arising: in the wake of a flapping foil have been observed in 
many experimental I Jones et q?Jll99S ; Godoy-Diana et a7ll2008llvon Ellenrieder fc PothosI 



f Note that the product of these two parameters gives the flapping amplitude based Strouhal 
number that is often used, i.e. St x Ad = fA/Uo = St a 
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Figure 2. Fluorescein dye visualization of the typical reverse Benard-von Karman vortex street 
that characterizes flapping based propulsion (top); and an asymmetric wake (bottom) that is 
produced by some flapping configurations even when the flap motion is symmetric. 



2008; Buchholz fe Smitsll2008h and numerical studies (jjones et a/Jl998tlLewin fe Haj-Hariril 



2003D even when the geometry of the problem and the flapping motion are symmetric. 



Other configurations of forced wakes have been shown to produce asymmetric flows that 
arise from the interaction of natural and forced vort ex shedding — see for i nstan ce the 
soap-film experiments with forced cylinder wakes in Couder fc Basdevant 1 19861 ). The 
deflection of the propulsive vortex street determines that the net force generated by the 
flapping motion is not aligned with the foil symmetry plane, or in other words that a 
mean lift force accompanies the production of thrust. A qualitative explanation of the 
process that determines the symmetry breaking can be given analyzing the structure of 
the wake: the deflection of the vortex street that signals the symmetry breaking results 
from the formation of a dipolar structure from each couple of counter-rotating vortices 
shed on each flapping period (see figure [2]). Above a certain threshold, the self-advection 
of the dipolar structure formed over one flapping period is strong enough to decouple 
from the subsequent vortex in the street and to generate a deflection of the mean flow. 
These asymmetric wakes occur in a region of the parameter space that overlaps the 



4 



R. Godoy-Diana, C. Marais, J.-L. Aider and J.E. Wesfreid 




^ measurement plane 

Figure 3. Schematic diagram of the experimental setup in the hydrodynamic tunnel. 



high-efficiency Strouhal number range used by flapping animals, which makes the precise 
definition of the symmetry breaking threshold potentially important for the design of ar- 
tificial flapping-based propulsors and their control. This is the goal of the present work, 
where we study in detail the spanwise vorticity field in the near wake of the flapping foil 
in order to characterize all the basic features of the reverse BvK vortex streets. Relying 
on a hypothesis of quasi-two-dimensionality of the flow in the near wake, we propose a 
predictive symmetry breaking criterion based on the phase velocity of the vortex street. 



2. Experimental setup 

The setup is the same described in lGodoy-Diana et al. ( 2008h and consists of a pitching 
foil placed in a hydrodynamic tunnel (see figure [3]). The foil chord c is 23mm and its span 
is 100mm which covers the whole height of the 100 x 150mm section of the tunnel. The 
foil profile is symmetric, opening at the leading edge as a semicircle of diameter D = 5mm 
which is also the maximum foil width. The pitching axis is driven by a stepper motor. The 
control parameters are the flow velocity in the tunnel Uo, the foil oscillation frequency 
/ and peak-to-peak amplitude A which let us define the Reynolds number, the Strouhal 
number and a dimensionless flapping amplitude as, respectively: 



Re = UD/v, St = fD/U , A D =A/D. (2.1) 

where v is the kinematic viscosity. In the strongly forced regimes produced by the flapping 
foil, the flapping frequency used to define St is equivalent to the main vortex shedding 
frequency. The boundary layer thickness on the tunnel walls at the position of the flap 
is of approximately 10mm for the present experiments. Measurements were performed 
using 2D Particle Image Velocimetry (PIV) on the horizontal mid-plane of the flap. PIV 
acquisition and post-processing was done using a LaVision system with an ImagerPro 
1600 x 1200 12-bit CCD camera recording pairs of images at ~ 15Hz and a two rod 
Nd:YAG (15mJ) pulsed laser. Laser sheet width was about 1mm in the whole 100mm 
x 80mm imaging region. The time lapse between the two frames (df) was set to 12ms. 
Additional p ost-processing and analysis were done using Matlab and the PIVMat Toolbox 
(|Moisvll2007h . 
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Figure 4. Three consecutive snapshots of the vorticity field in the wake of the flapping foil at 
Re — 255 and Ad = 1-42 for two different flapping frequencies. Left column: St = 0.15, reverse 
BvK wake; and right column: St = 0.27, asymmetric wake. 
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Figure 5. Circulation (a,e), radius (b,f), x-position (c,g) and y-position (d,h) of two consecutive 
(and counter-rotating) vortices in the wake of the flapping foil as a function of time. The top 
figures (a-d) correspond to the reverse BvK wake and the bottom figures (e-h) to the asymmetric 
wake shown in figured] (left and right columns, respectively). 



3. Observations 

3.1. The vorticity field 

The lo z vorticity field is calculated from the u x and u y fields obtained from PIV using 
2nd-order centered differences. Two sequences of vorticity fields are shown in figured! the 
left column corresponds to a reverse BvK street and the right column to an asymmetric 
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Figure 6. Definition of the asymmetry angle 6 from the deflected mean jet flow (left) and 9 as 
a function of the Strouhal number for Ad = 1.4 and 2.1 (right). The circles shown in the left 
plot are the maxima of cross-stream profiles of the average horizontal velocity field. A linear fit 
over these points in the near wake gives the dotted line that defines the angle theta. The error 
bars in the right plot represent a constant ±1.5 degrees that was the maximum deflection angle 
measured with the previous method for a wake in the symmetric reverse Benard-von Karman 
configuration. 



wake. The Reynolds number and napping amplitude are the same in both cases and 
only the Strouhal number has been increased from 0.15 for the reverse BvK to 0.27 for 
the asymmetric wake. The positions of two consecutive (and counter-rotating) vortices 
are followed for each experiment using a search of local maxima (X max (t),Y max (t)) and 
minima (X min (t), Y min (t)) in the lu z field. These are the coordinates of the center of the 
rectangles in figure [H 

3.2. The circulation T 

An area of integration that encompasses each vortex needs to be defined in order to 
calculate its characteristic features such as circulation and size. We use gaussian fits 
cxp(— xf /of ) along the vertical and horizontal axes centered on the positions of the 
maxima and minima of vorticity and define the sizes of the vortex along each direction 
x and y as 2<Tj . These are the horizontal and vertical sizes of the rectangles in the time 
sequences of figure [H The choice of rectangular integration contours (instead of the 
elliptical ones that would have better followed the vortex shape) was kept in order to 
avoid introducing error from interpolation of the PIV data. Since the vortex cores are 
nearly circular we define a single vortex radius a as the mean of the sizes calculated 
along the two principal axes. The circulation T can be then calculated either from a 
line integral of the velocity field on each contour or a surface integral of the vorticity 
field over each rectangle. Although in theory the two definitions of T are equivalent, 
the two calculations from the experimental fields are slightly different mainly because, 
vortices in the wake being not too far from each other, the contour for the line integral 
calculation sometimes passes through a neighboring opposite-signed vortex (see figure 
SI right column) and gives thus a spurious contribution to the total circulation. The 
circulation plotted in figure [5] is the mean value between the two different methods of 
calculation, the difference giving the error bars. When the integration regions for the 
two counter- rotating vortices do not overlap, as is the case in the left column of figure [H 
the difference in the two calculations is small (see figure 03a). On the contrary, the error 
bars are larger when the two rectangles overlap, as in the case of the asymmetric wake 
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Figure 7. (a) Circulation To and (b) distance between two consecutive vortices £o as a function 
of St and Ad- Circles, triangles and squares correspond, respectively, to Ad = 0.7, 1.4 and 2.1. 



depicted in the right column of figureHl This can be seen in figure[5je, where the spurious 
effect due to having the neighboring vortex partially overlapping the integration region 
appears as a decreasing trend (in absolute value) in the time evolution of T. We also plot 
in figure [5] the time evolution of the vortices radii and positions. It can be seen in figure 
[5jh that the deflection of the wake is correctly captured by the time evolution of the 
y-coordinate. The time axes in all these plots are nondimensionalized by the advection 
time scale r = D/Uq. 



3.3. The asymmetric wake 

The domain of existence of the reverse BvK vortex street is bounded on the upper-right 
zone of the [St, Ajj) phase space by a transition to an asymmetric regime (figure [1]). 
In order to characterize this transition we define an angle of asymmetry 9 using the 
direction of the jet observed in the mean velocity field (see figure E]). Of the three series 
with different flapping amplitudes studied here (see figure [7| in the next section) , the 
smallest one (Ad = 0.7) does not produce an asymmetric wake so only the cases of 
Ap = 1.4 and 2.1 are reported in figure [6l These measurements show that the transition 
is rather abrupt and support the idea of the existence of a symmetry breaking threshold. 



4. Discussion 

4.1. Back to the (St, Ad) phase space 

Thinking in a (St, Ad) phase space perspective, we use the 'initial' values of the cir- 
culation and the positions of the maxima and minima of oj z , defined as the values 
measured when the positive vortex crosses a reference frame fixed at X Ic i = 2.5D, in 
order to compare the vortices produced by different flapping configurations. We thus 
plot in figure [7| the initial value of the circulation Tq and the distance between the 
two consecutive counter-rotating vortices £o , as a function of St for three series with 
different values of Ad- The definition of £ in terms of the positions of two counter- 
rotating vortices shown in figure [U is, as shown graphically in figure [9] (left), £ = 
\Z(X max - X mm ) 2 + (Y max - Y min ) 2 . The main observation is that, for a given ampli- 
tude, the circulation increases and the distance between consecutive vortices decreases 
with increasing flapping frequency. Although this is not surprising, it is an important 
point because it lets us ascertain that a threshold curve determining the symmetry 
breaking of the reverse BvK wake can be traced in a (ro,£o) space and also that a 
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Figure 8. (a) Phase velocity [/phase of the vortex street normalized with respect to the 
free-stream Uo; and (b) translation speed of the initial dipoles [/dipole calculated using Eq. 
I4.ll Both quantities are plotted as a function of St and Ad- Circles, triangles and squares cor- 
respond, respectively, to Ad = 0.7, 1.4 and 2.1. Asymmetric wakes are marked with a dark 
dot. 



model containing the basic physics of the problem could be tested using experimental 
measurements. 

4.2. Phase velocity 

The advection velocity of the vortices with respect to the free-stream can be calculated 
from the re-position measurements Xt(t) in figuresEJc andODg. We define this phase speed 
[/phase as the slope dXi / dt evaluated at the initial reference time Q. In the case of the 
Benard-von Karm an vortex street past a bluff body [/phase < Uq (see e.g. lBearmanlll967l : 
IWilliamson 1989f), which is related to the velocity deficit in the wake of the obstacle. The 



phase velocity for the present experiments is plotted in figure [8] (left) . The propulsive 
nature of the reverse BvK wake appears clearly in the fact that [/phase > Uq. 

4.3. The dipole model 

The physical mechanism giving rise to a deflected wake is based on the formation of a 
dipolar structure on each flapping period ( Godoy-Diana et all 2008|). a feat ure that has 



also been observed in forced wakes in soap films ( Couder fc Basdevant 198fJ ) , supporting 



the idea of a mainly 2D phenomenon. The initial condition sets the choice of the side 
where the asymmetry develops: the first dipole that is formed entrains fluid behind it, 
deflecting the mean flow in the wake and forcing the subsequent dipoles to follow the same 
path. This initial perturbation exists for all the (St, Ajj) parameter space, however, after 
a few periods of symmetric flapping motion, only in a region above a certain threshold 
(see figure[l]) the wakes remain asymmetric. The two quantities (Fo,£o) plotted in figure 
[7] can be used to give a measure of the strength of the dipolar structures that are formed 
for each flapping configuration. The simplest model that contains T and £ is to consider 
a dipole made of two point vortices of circulations ±F separated by a distance £. In this 
case, the translation speed of the dipole ( determined by the effect of each vortex over the 



other) is given by (see e.g. ISaffma: 



lipole ( d 
nl ll992l l 



f Vortices slightly accelerate in the initial part of the wake so that the wavelength measured 
in a snapshot of the vorticity field over the vort ices nearest to the flap is shorter than the 
wavelength measured farther in the wake (see also lBearrnanj[l967)- Here we deliberately define 
[/phase as the advection speed at the initial reference time because we are interested in the near 
wake mechanism that triggers the symmetry breaking. 
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Figure 9. Left: Definition of the angle a used in eauation l4.2l The line at angle a with the hor- 
izontal is at right angles with the segment of length £ that joins the two consecutive counter-ro- 
tating vortices. Right: Effective phase velocity U p * normalized with respect to the free-stream 
Uo plotted as a function of St and An- Circles, triangles and squares correspond, respectively, 
to Ad = 0.7, 1.4 and 2.1. Asymmetric wakes are marked with a dark dot. The points above the 
dashed axis verify the proposed symmetry breaking condition U* > 0. 



t/dipole - 2< 



(4.1) 



The values of [/dipole calculated using the data from figure [7| are shown in figure [8] 
(right). As a direct consequence of the behavior of Tq and £o in figure[7J for each flapping 
amplitude, the self-induced speed of the dipolar structure increases with the Strouhal 
number. It is remarkable that [/dipole can reach values up to almost 50% of Uq. The two 
plots in figure [8] show that a correlation exists between [/dipole and [/phase- This can be 
easily understood considering that the flow field produced by each vortex, and hence 
E^dipoie, contributes to the overall advection velocity in the vortex street. 

Thinking about the definition of a criterion for the symmetry breaking, it seems evident 
from figure[8]that no single threshold can be established for [/dipole or [/phase that accounts 
for all observations of asymmetric wakes (we have marked these in figure [8] with black 
dots). Because we have observed that the symmetry breaking is related to the ability 
of a given dipolar structure to escape from the "symmetrizing" effect of the subsequent 
vortices in the wake, it is useful to find a quantitative measurement of the velocity at 
which vortices are moving in the direction defined by the dipole. In order to do so we 
define the effective phase velocity, in terms of the angle a between the horizontal (the 
streamwise direction) and the direction of [/dipole (see figure [9]), as: 



U* = ([/phase - Uq) COS a - [/dipole ■ (4.2) 

Recalling that the measured phase velocity [/phase results from the superposition of the 
free-stream velocity and all the velocities induced by the vortices in the wake, the effective 
phase velocity defined in equation 221 is actually the component along the dipole direction 
of the velocity induced by the vortex street excepting the contribution of the dipole being 
considered, and not including the free-stream velocity. Within the limits imposed by the 
experimental uncertainty, it can be seen in figure [9] that a fair prediction of the observed 
behavior can be given saying that the reverse BvK vortex street will be prone to breaking 
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symmetry if the effective phase velocity U* is positive, which is when the vortex street 
downstream of the dipole being considered enhances the dipole velocity. 



5. Concluding remarks 

Reverse Benard-von Karman (BvK) vortex streets are a fundamental feature of fish- 
like swimming and forward flapping flight and their symmetry properties are intimately 
related to the cycle of thrust and lift production. Although the deflection of these propul- 
sive wakes has been observed and characterized in various flapping- foil experiments 
[ Jones et a/1ll998tlGodov-Diana e/:"alll2008tl7on Ellenrieder fc Pothosll2008h and numer- 
ical simulations ijLewin fc Haj-Haririll2003f h the model proposed in the present paper is, 
to our knowledge, the first attempt to produce a quantitative threshold prediction based 
on the observed physical mechanism underlying the symmetry breaking. In spite of the 
strongly idealized dipolar model used in Eq. l4.ll as a fundamental element of the model, 
the predicted threshold accounts reasonably well for the experimental reality and should 
be useful to rationalize the similar observations reported in the literature. 

We gratefully acknowledge Maurice Rossi for fruitful discussions on the dipole model. 
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